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Abstract 

In this article, we consider linearly stable elliptic fixed points (equi- 
librium) for a symplectic vector field and we prove generic results of 
super-exponential stability for nearby solutions. We will focus on the 
neighbourhood of elliptic fixed points but the case of linearly stable 
isotropic reducible invariant tori in an Hamiltonian system should be 
similar. 

More specifically, Morbidelli and Giorgilli have proved a result of 
stability over super-exponentially long times if one consider an analytic 
lagrangian torus, invariant for an analytic hamiltonian system, with 
a diophantine translation vector which admit a sign definite torsion. 
Then, the solutions of the system move very little over times which are 
super-exponentially long with respect to the inverse of the distance to 
the invariant torus. 

The proof is in two steps: first the construction of a Birkhoff nor- 
mal form at a high order, then the application of Nekhoroshev theory. 
Bounemoura has shown that the second step of this construction re- 
mains valid if the Birkhoff normal form linked to the invariant torus 
or the elliptic fixed point belongs to a generic set among the formal 
series. This is not sufficient to prove this kind of super-exponential sta- 
bility results in a general setting. We should also establish that most 
strongly non resonant elliptic fixed point or invariant torus in a Hamil- 
tonian system admit a Birkhoff normal form in the set introduced by 
Bounemoura. We show here that this property is satisfied generically 
in the sense of the measure (prevalence) through infinite-dimensional 
probe spaces (that is an infinite number of parameter chosen at ran- 
dom) with methods similar to those developed in a paper of Gorodetski, 
Kaloshin and Hunt in another setting. 
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1 Introduction 



We are interested in the stability properties, in the sense of Lyapounov, of 
linearly stable elliptic fixed points (equilibrium) denoted on a symplectic 
manifold (M,Cl) for a symplectic vector field X (i.e. : ix^ is closed), that 
is X(x*) = 0. 

As the problem is local, it is enough to consider a Hamiltonian H, defined 
and analytic on an open neighbourhood of £ 1" x R n equipped with the 
standard symplectic form, having the origin as a fixed point. Up to an 
irrelevant additive constant, expanding the Hamiltonian as a power series 
at the origin we can write 

H(z) = H 2 (z) + V(z) 

where z = (x, y) € W 1 x ]R n is sufficiently close to the origin in W l x W\ H 2 
is the Hessian of H at and V(z) = 0(|z| 3 ). 

Recall that the fixed point is said to be elliptic if the spectrum of the 
linearized system is purely imaginary. Due to the symplectic character of the 
equations, these are the only linearly stable fixed points and the spectrum 
has the form {zkiui, . . . , ±ioj n } for some vector oj = (ui, . . . ,oj n ) which is 
called the normal (or characteristic) frequency. Now we assume that the 
components of oj are all distinct so that we can make a symplectic linear 
change of variables that diagonalize the quadratic part, hence 

H(z) =lo.1 + V{z) 

where 1 = l(z) are the "formal actions" , that is 

1 n 
X{z) = (Zi, . . . ,l n ) = -{x\ + yl . . . , x 2 n + yl) e R n and uj.X = Y^Pr 

i=i 

If all the components of oj have all the same sign, then H is a Lyapounov 
function and the fixed point is stable. But in a general setting, one has to 
take into account the higher order terms V(z) = 0(jz| 3 ). 

Given a solution z(t), if we denote Z(t) = Z(z(t)), then ||Z(t)|| is es- 
sentially the distance of z(t) to the origin and \\Z(t) — Z(0)|| measures the 
deviation of z(t) from its original position. 

Without loss of generality, we can assume that we consider a Hamiltonian 
defined on the ball Br of radius R < 1 in R 2n . Then by analyticity, we 
extend the resulting Hamiltonian as a holomorphic function on some closed 
complex ball D s of radius s < 1 in C 2n . If we define A s as the space of 
holomorphic Hamiltonians on D s which are real valued for real arguments, 
and ||.|| s its usual supremum norm, then we are led to consider 

(H(z)=cu.l+f(z) 
[H e A s ; f(z) = 0(z*). 
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To obtain results of stability, one of the main tool is to construct normal 
forms via classical averagings, and in this case these are the so-called Birkhoff 
normal forms. 

For an integer m > 1, assuming w is non-resonant up to order 2m, that is 

n 

fc.w/O, k E Z n , < p||i =^2\kj\ < 2m 

i=i 

then there exists an analytic symplectic transformation & m close to identity 
such that H o & m is in Birkhoff normal form up to order 2m, that is 

m 

H o $ m (z) = h m {l) + f m (z) where h m {X) = £> (fc) (^) 

k=l 

with homogenous polynomial of degree k in the X variables and the 
remainder f m of order z 2m (see [Bir66] or [Dou88j, [HZ94 j for a more recent 
exposition). The polynomials (B ^ fc ') fcgN » are uniquely defined (even though 
the transformations are not) and are usually called the Birkhoff invariants. 
The transformed Hamiltonian is therefore the sum of an integrable part h m , 
for which the origin is trivially stable in the sense that X(t) is constant for 
all times, and a much smaller perturbation f m for a small enough neigh- 
bourhood of the origin. Moreover, if u is non-resonant up to any order, we 
can even define a formal symplectic transformation and a formal power 
series /ioo = ^# (fc) such that 
fe>i 

Ho$ 00 (z) = h 00 (X). (1) 

However, for a certain topology on the coefficients, the transformation 
^oo is generically divergent as it was proved by Siegel |Sie41j and the con- 
vergence properties of the series are even more subtle. Actually, Perez- 
Marco ( |PM03| ) has showed that if h oo is divergent for some H G <A S , then 
the complete Birkhoff normal form hoo is divergent for a "typical" H G A s - 
But nevertheless, we cannot ensure that is generically convergent or 
divergent. 

The first kind of stability result is given by an application of KAM theory. 
Assume that uj is non-resonant up to order 4 so that our system reduce to 

H(z) = h 2 (z) + f 2 (z) = u.X + B^(X) + f 2 (z). 

where is a quadratic form. 

Under the so-called Kolmogorov non-degeneracy condition (or the twist 
condition in the context of an elliptic fixed point) : the quadratic form B^ 
is non- degenerate, KAM theory (see [AKN97J and [DG96] for the specific 
case of an elliptic fixed point) ensures stability for all times for most of the 
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initial conditions. That is, for a large set of solution z(t) : if |Z(0)| < p 
then the variations \T{t) — X(0)| are of order p for all times, hence we have 
Lyapounov stability. More precisely, KAM theory yields the existence of 
invariant Lagrangian tori which form a set of positive measure with a density 
going to one at the origin. However because of the dimensions, it is only for 
n = 2 that we can deduce perpetual stability results. Actually, for n = 2 the 
invariant tori are 2-dimensional and divide the 3-dimensional energy level, 
therefore the solutions of the perturbed system are global and bounded over 
infinite times. On the other hand, an arbitrary large drift of the orbits can 
still occur for n > 3 in the complement of the invariant tori which is a 
dense set with a complicated topology. In fact it is believed that, except 
for two degrees of freedom systems (and of course one degree), "generic" 
elliptic fixed points are unstable (see |DLC83| and [Dou88j for examples and 
[KMV04] for an announcement in a generic setting). 

Thus, for n > 2, KAM theory provides results of stability only in the 
sense of measure and theorems of stability which are valid for an open set 
of initial condition can only be proved over finite times. 

Definition 1.1. The origin is stable over finite times if \I(0)\ < p small 
enough implies 

\X{t) -X(0)| = R{p) for \t\ < T{p) 

where R(p) is of order p and T(p) is at least of order p~ 1 (otherwise we 
have a trivial bound). 

We have polynomial stability ifT(p) is at least of order p~ m for some 
fixed integer m G N* . 

We have exponential stability ifT(p) is at least of order exp (Cp~ a ) 
for some fixed constants C > and a > 0. 

We have superexponential stability if exp (Cp~ a ) is negligible with 
respect to T(p) for any constants C > and a > 0. 

Using the mean value theorem, a polynomial time of stability can be 
easily obtained once a Birkhoff normal form up to a finite order has been 
built. The first result of this kind around an elliptic fixed point was obtained 
by Littlewood in 1959 f [Lit59j ). 

To obtain results of exponential stability, there is basically two distinct 
methods. 

In the first approach, we call Birkhoff type estimates, one assumes that 
the frequency uj is a (^,r)-Diophantine vector for some positive constants 7 
and r, hence u) G 7jT - where : 

fi 7 , r = La e M n such that |Jfe.w| > -pyp- for all k G Z"\{(0, ... ,0}| (2) 

We recall that the measure of the complementary set of 7 T is of order 
0(7) for r > n - 1. 
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In particular u> is completely non-resonant, hence one can perform any 
finite number m of Birkhoff normalizations, and since we have a control on 
the small divisors, one can estimate the size of the remainder f m . 

Under these assumptions, one can prove ([B G"86] , [GDF + 89] ) that for a 
real-analytic system, the "action variables" I become quasi integrals over 
exponentially long times, more specifically : 

Theorem 1.2. Consider an elliptic fixed point for a real analytic Hamilto- 
nian system of type (A) . 

There exists positive constants C\, C 2 , C3, C4 which depend only of the 
analyticity width s, the size ||-ff|| s of the Hamiltonian H, the number of 
degree of freedom n and t the exponent in the Diophantine condition such 
that if p < C17, an arbitrary solution (x(t),y(t)) with an initial action 
11^(0)11 < p is defined at least over an exponentially long time and satisfies : 

\\l(t) - J(0)|| < C 2P if \t\ < C 3 exp (c 4 (£) 1+T ^j . (3) 

The proof rely on the existence of a normalizing transformation up to 
an exponentially large order m comparable with the inverse of the distance 
to the origin. This is possible since we have lower bounds on the small 
denominators and the growth of the coefficients in the normalizing expan- 
sion is reduced to a combinatorial problem (the "great multipliers" along 
Poincare terminology) . Actually, under the assumptions of linear diophan- 
tine frequencies and the real analyticity of considered Hamiltonian, there 
exists a decreasing sequence of radii s m ~ j/m 1+T for m € N* such that 
the normalizing transformation and the Birkhoff normal up to order m 
can be build on the ball of radius s m with a remainder of size compara- 
ble with s™.(m\) 1+T ~ (m\/m m ) 1+T which is dominated by exp(— (1 + r)m) 



-ex P (-C(^)-). 



Finally, since the averaged Hamiltonian is integrable, the speed of drift 
of the action variables is very slow. 

More specifically, there exists a sequence of annuli of decreasing radii 
(s m )meN* such that : 



Wi < I |Z(0) 1 1 < s m ^\ \T{t) - 1(0) \\<C 2 s m if |i|<C : 




In order to obtain uniform estimates with respect to the initial radius 
p, as in our theorem 1.2, it only remains to slightly lower the coefficients 
C2,Cs,C4. Treschev continuous averaging method ([Tr02]) would allow to 
obtain directly a result with respect to the initial radius p and to avoid 
this sequence (s m ) mS N* but we will stay in the setting of Giorgilli et al 
( |GDF+89j . see also |DG96j ). 
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The second approach, we call Nekhoroshev type estimates, is fundamen- 
tally different. One just assume that cj is non-resonant up to order 4 so that 
the Hamiltonian reduces as for KAM theory to 

H(z) = h 2 (z) + f 2 (z) = uj.I + B^{X) + f 2 (z). 

where B^ is a quadratic form. But this time we require that the non-linear 
part h 2 (z) = uj.X + B^\X) is convex in actions X, which is equivalent to B^ 
being sign definite. Under those assumptions, it was predicted by Lochak 
( |Loc92| ) and proved independently by Niederman ( |Nie98| ) and Benettin, 
Fasso and Guzzo ( |BFG98j ) that exponential stability holds. Their proofs 
are based on the implementation of Nekhoroshev estimates in cartesian co- 
ordinates, but they are radically different : the first one uses Lochak method 
of periodic averagings and simultaneous Diophantine approximations while 
the second is based on Nekhoroshev original mechanism. The proof of Nie- 
derman was later clarified by Poschel ( |Pos99j ). 

Actually, Birkhoff estimates rely on the existence of an accurate inte- 
grable normal form while Nekhoroshev theory (see [Nek77j or |Niellj for a 
survey) gives exponential stability for all the solutions even if resonances of 
low order are present. This problem is overcome by a suitable transversality 
condition called steepness on the geometry of the integrable system. Steep- 
ness is a generic condition infinitely differentiable functions which is implied 
by (quasi)-convexity. 

A remarkable result of Morbidelli and Giorgilli QMG95J) clearly shows 
that the two previous results which come respectively from Hamiltonian 
perturbation theory and Nekhoroshev's theorem are independent and can 
be superimposed. 



Theorem 1.3. Consider an elliptic fixed point for a real analytic Hamilto- 
nian system of type (A) with a Diophantine normal frequency uj and a sign 
definite torsion, that is : B^ is a sign definite quadratic form. 

There exists positive constants C\, C 2 , C3, C4 which depend only of 
the analyticity width s, the size \\H\\ S of the Hamiltonian H, the number 
of degree of freedom n and r the exponent in the Diophantine condition 
such that if p < C17, an arbitrary solution (x(t),y(t)) with an initial ac- 
tion ||X(0)|| < p is defined at least over a superexponentially long time and 
satisfies : 

\\X(t) - J(0)|| < C 2P if |i| < C 3 exp ^exp (c 4 {^j (4) 

The proof starts with a Birkhoff's normal form with an exponentially 
small remainder like in the theorem 1.2. Then, Nekhoroshev's theory is 
applied with a perturbation which is already exponentially small, hence we 
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obtain a superexponential time of stability. More specifically, we need the 
(quasi-)convexity of the Birkhoff normal forms at all order. Since (quasi-) 
convexity is an open condition for the —topology, the convexity of the 
second Birkhoff invariant implies the same property for h m (Vm € N*). 

One of the main problem to extend the later kind of results in a generic 
setting is that we need a property on the Birkhoff invariants at all order 
which is not an open condition. 

In this direction, Bounemoura ([Boulla]) has proved that these reason- 
ings remain valid in the case of a non resonant elliptic fixed point under a 
generic condition on the associated Birkhoff normal form /i^, more specifi- 
cally : 

Theorem 1.4. Consider an elliptic fixed point for a real analytic Hamilto- 
nian system of type {A) with a non-resonant normal frequency oo. 

If its Birkhoff normal form hoo belongs to a generic set in the space of 
the formal power expansions M[[X]] (genericity is in the sense of prevalence, 
see the next section) the estimates of superexponential stability given in the 
previous theorem remain valid. 



2 Setting and statement of our results : 

Bounemoura gave only a partial answer to the problem of superexponential 
stability of elliptic fixed position. Actually, we cannot ensure that generically 
the Birkhoff normal form associated to a non resonant elliptic fixed point 
belongs to the set introduced by Bounemoura even if it is a big set in the 
sense of measure. This problem is not straightforward since we have very few 
general information on the Birkhoff normal forms. Actually, in the context 
of KAM theory, we are only aware of one result of Eliasson, Fayad and 
Krikorian (so far unpublished) which involves all the Birkhoff invariants and 
ensure the existence of invariant tori which accumulate an elliptic fixed point 
in a Hamiltonian system with a Diophantine linear spectrum but without 
any torsion condition. 

Here, we provide a complete answer to our problem of stability. 

We first specify our space of Hamiltonians 



H(x, y) = ijj.X + f(x, y) where 



H € A s , 

u £K" and f(x,y) = 3 (x,y). 



we recall that A s is the space of holomorphic Hamiltonians on D s C C 2n the 
ball of radius s around the origin which are real valued for real arguments 
equipped with the supremum norm ||.|| s on D s , hence it is a Banach space. 



Without loss of generality, we can assume that s < 1. 
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For an homogenous polynomial P(I) = PkZ k , as in [GHK06] we 

\l\=k,lGN 2n 

consider the following Euclidean norm (also called Bombieri norm) 

\i\=k k 

with the multinomial coefficient C\. 

This norm is invariant by isometry in C 2n (see [GHK06J). 

There exists a Borelian measure canonically associated to the norm | \P\ \k 
on the space Vk of homogenous polynomials of degree k. 

We consider the Hilbert brick 

HB N = J h(l) = h k (X) with \\h k \\ k < 1 (Vfc) 1 
I ken* J 

it is a compact set for the product topology when the space of polynomial 
is identified as the infinite product J^J Vk- 

fceN* 

This compact can be continuously embedded in our Banach space of 
Hamiltonian 7i s with the assumption that s < 1, hence it is a compact set 
in (H s , \\.\\s). 

If we denote by pk the Borelian probability measure on the unit ball of 
V k for all k G N*, the embedding of HB N in (U s , \\.\\ s ) can be equipped 
with the image of the product measure [i = /ifc which define a compactly 

kern* 

supported Borelian probability measure in (T-L s , ||.|| s ). 

This setting was previously considered by Gorodetski, Hunt and Kaloshin 
([GHK06]). 

Theorem 2.1. (Main theorem) For an arbitrary H E Us, the modified 
Hamiltonian H + h is super- exponentially stable for almost all h G HB N 
with respect to the probability [i. 

More precisely, there exist positive constants a, C\, C2, C3 such that for 
a small enough p, every solution of H + h with initial condition \I(0) \ < p 
satisfies 

\X{t) -X(0)| < dp for \t\ < exp [c 2 p- c ^ ln ^l") . 

In other words, the set of super exponentially stable Hamiltonians is 
prevalent in our Banach space since p defines a compactly supported Bore- 
lian measure which is transverse to the set of non super-exponentially stable 
Hamiltonian (see |()Y05j and [HKlOp . 

We do not see obstructions to obtain similar results for Lagrangian tori 
(see the remarks at the end of this section). 
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3 Scheme of the proof : 



Either in this paper or in Bounemoura work ( [Boulla] ). the proofs rely 
crucially on the extension of Nekhoroshev theory obtained in |Nie07] for 
a much wider class of unperturbed Hamiltonian than the steep integrable 
Hamiltonian considered originally in [Nek77]. 

For instance, it can be proved ( [Boulla] ) that almost all quadratic form 
in action variables fall into this class of Diophantine Morse functions while 
only sign definite quadratic forms are steep according to Nekhoroshev defi- 
nition. 

Actually, the definition of Diophantine Morse functions involves count- 
able quantitative transversality conditions which are stated in adapted co- 
ordinates. It is inspired on one hand by the steepness condition introduced 
by Nekhoroshev ( [N ek77j ) and on the other hand by the quantitative Morse- 
Sard theory of Yomdin (|Yom83j and [YC04]) where we consider "nearly- 
critical" points which are quantitatively non degenerate. 

The application of Yomdin quantitative Morse Sard theory allows to 
prove ( [Nie07| ) that for an arbitrary real analytic integrable Hamiltonian h 
defined in a neighbourhood of the the origin in the action space M n , then 
for almost any parameter u £ M. n , the modified integrable Hamiltonian 
huj(x) = h(I) —uj.I is exponentially stable (i.e. : Nekhoroshev estimates are 
valid) . 

In the case of an elliptic equilibrium, we have to deal with cartesian co- 
ordinates instead of action-angles ones and we cannot use the latter since 
they become singular at the origin, this yields specific problems detailed 
in [Loc92j and overcomes in [BFG98J, [Nic98] for an elliptic point with a 
sign definite torsion. In [BN12], Nekhoroshev estimates are proved in a 
generic setting by means of successive periodic averagings and this method 
can be implemented in cartesian coordinates. These observations allowed 
Bounemoura [BoullaJ to prove that a Morse Diophantine condition is satis- 
fied at all order on a full Lebesgue measure set of Birkhoff invariants which 
is strong enough to prove super exponential results of stability over times of 
order exp(exp(C/r a )). 

As it will be specified in the sequel, the set considered by Bounemoura is 
not big enough to ensure that generically the Birkhoff normal form of a non 
resonant elliptic fixed point falls into this good set. We will have to consider 
more refined estimates on the set of bad parameter oj € W 1 according to the 
following proposition which is a direct extension of ( |Nie07| . theorem 3.2.6) 
or f |BJN12j . theorem 2.2) : 

Proposition 3.1. Consider an integrable Hamiltonian 
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holomorphic on the ball of radius s around the origin in C . 

There exists an open set Q C B p where B p is the ball of radius p around 
the origin in W 1 such that 

Vol(C^) < C p e a 

where C p > depends only of n, s, p, ||H|| S and a > depends only of n. 

For uj G B p \Ce P \ Nekhoroshev estimates of stability can be proved for 
the modified integrable Hamiltonian H^x, y) = h w (I) = h(Z) — uj.Z with an 
holomorphic perturbation f(x,y) of size \\f\\ s < £■ 

Remark 3.2. Note also that h does not need to be polynomial, we only 
require that the Hamiltonian is smooth enough. 

Now, following the strategy of Morbidelli and Giorgilli QMG95J), the 
estimates of superexponential stability given by Bounemoura rely on the 
existence for all m G N* of a Birkhoff normalization $ m specified in theorem 
1.2. which integrate H up to order 2m on the ball of radius s m > around 
the origin in C n , hence H o $> m (z) = h m (Z) + f m (z) where h m G M m [X] and 
the remainder f m is of order z 2m . Moreover, the size of the perturbation 
ll/m||s m \ s compar able with exp (-(1 + r)m) ||/i|| Sl and ||/i|| Sl ~ si (see 
[GDF+89] . [DG96j ^. 

Then, in order to obtain superexponential estimates of stability, we 
should be able to apply Nekhoroshev theory (with the proof given in [BN12J) 
on the Birkhoff normal forms (/i m ) mgN * at all order m with a perturbation 
f m which decrease geometrically. 

Going back to the genericity problem, according to the previous propo- 
sition for an arbitrary m G N*, Nekhoroshev estimates can be proved on the 
ball of radius s m > for the perturbation f m of h m (Z) + uj.Z if to G B P \C$ 

where Cm^ is an open set which contain the critical parameters in the 
ball B p C K n with a volume which decrease geometrically with m hence 

vol ( U C ^ ] ) < E Vo1 ( c ^ p) ) - Vo1 { c i ] ) < °°- 

VmeN* / mGN* 

In our infinite-dimensional dimensional setting, we will say that a formal 
expansion hoc G is superexponentially stable if there exists a sequence 

of radii (s m ) mGN , such that Nekhoroshev estimates can be proved for all 
mGN* with the integrable part given by the truncation Proj Rm r£i (/loo) at 
order m and a perturbation f m of size ||/ m || Sm ~ si/e <yl+T ^ m . 

We note that we have Proj^^/ioo + t u ) = Proj ffim[X ] (hoc) + L where 
£u(X) = uj.X for all cu G K n hence, for an arbitrary formal expansion /i^ in 
the modified expansion hoo(I) + ui.Z is superexponentially stable if 
ueB p \\J d p) . 

meN* 
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Prom these estimates Bounemoura ( [Boulla] ) proves that, for an arbi- 
trary formal expansion € superexponential estimates of stability 
are valid for h QO (Z) - w.l on a ball of small enough radius e around the 
origin if co € B p \C { £ p) where C { £ p) = (J C$ with Vol of the order of e a 

meN* 

where a given in the previous proposition depends only of n. 

This implies the generic property stated in Theorem 1.4. 

Especially, Bounemoura QBoullaJ) shows the existence of a foliation 
of the space of formal series M[[X]] in n-dimensional probe spaces directed 
by the space of linear forms where the properties required to prove super- 
exponential estimates are satisfied at almost every point on the leaves for 
the n-dimensional Lebesgue measure. The issue is that when we want to 
come back in the initial (non- normalized) coordinates, a change of the lin- 
ear frequencies implies a change of all the Birkhoff invariants and all the 
normalizing transformations. Hence it is impossible to recover a foliation 
of the initial space of Hamiltonian H s by n-parameter families where su- 
perexponential stability can be proved for almost all parameters with the 
n-dimensional Lebesgue measure. We must use more global estimates with 
an infinite number of parameters and the important observations in this 
direction are the following. 

To state our results we have considered the embedding of HB^ in the 
Banach space 7i s . Here, we first consider the Hilbert brick HB 1 ^ embedded 
in the space of formal expansion identified with the product space 

n^*^ equipped with the product topology (we recall that Vk is the space 
of homogenous polynomial of degree k € N*). 

Actually, despite we are in an infinite dimensional space, the central es- 
timates concern polynomial integrable Hamiltonian (the truncated Birkhoff 
normal forms) and now we revisit the study of Bounemoura for polynomial 
Hamiltonian. For all m € N*, we identify the space M m [X] of polynomial of 
degree bounded by m with the product J~| Vk of the spaces (Vk)keN* of 

l<fc<m 

homogenous polynomial of degree k < m equipped with the product topol- 
ogy. Moreover, for all m € N*, we denote the projection Proj R rjfi (i?i? N ) 

by HB- m C R m [X] equipped with the product measure fj,< m = 

l<k<m 

We note that, for all m € N*, we have \i = /x< m x j^J fik and especially 

k>m+l 

fj, = fj,i x JJ/Xfc where fj,\ is the classical Euclidean volume up to a scaling 

k>2 

(to define a probability) on the unit ball in the space of linear forms. 

As in the previous reasonings, for all m £ N*, we consider an open set 
Xi m) in HB^ m C R m [X] given by Yomdin Morse- Sard theory which include 
the expansions in HB- m (considered as integrable Hamiltonians) where 
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Nekhoroshev theory cannot be applied to prove superexponential results of 
stability on a ball of radius e. Using Bounemoura construction, we can 
prove that intersects any affine subspace in lR m [X] directed by the 

vectorial space of linear forms along a set of volume dominated by (e/^ m )° 
with a > 0, given in proposition 3.1, which depends only of the number of 
degree of freedom n and a constant £ > 1. 

m 

With the fact that pL< m = fJ>i x JJ^fc where n\ is the Euclidean volume 

k=2 

771 

and JJ/^fc is a probability measure, the application of Fubini theorem implies 

k=2 

that the measure n< m (^X^ J remains is of the order of 

Now we come back to our initial functional space 7i s , for an arbitrary 
Hamiltonian H G H s we consider the application : 

BMT m : M m — ► R m [X] with R m [X] ~ H V k 

l<k<m 

(Pl,...,P m ) l > (Ql,...,Q m ) 

where Qk is the Birkhoff invariant of order k £ {1, . . . , m} of H + P for 

m 

P = and .A4 m is an open set in M m [JC] where the Birkhoff normal 

fc=i 

form at order m can be build. 
We have : 

Proposition 3.3. For k = 1, the transformation is Q\ = P\ + B^'(H) 
where (H) is the first Birkhoff invariant of H hence it is a translation. 
For k G {2, . . . , m} : 

Qk = Pk+ Fk(Pi, • • • , Pfc-l) 

where Tk is a differentiable function over M. m , hence the mapping BNT m 
is one-to-one with a Jacobian equal to 1. 

A similar property on the Newton interpolation polynomials was used 
by Gorodetski-Hunt-Kaloshin ( |GHK06j ). 

Despite our transformation has a unit Jacobian, it does not preserve the 
measure ^< m since it does not leave invariant the projection HB- m of the 
Hilbert Brick HB n . 

Actually, for an analytic Hamiltonian, the Birkhoff normal form should 
admit a Gevrey regularity and not remain analytic (we recall that we have 
no results on the generic divergence or convergence of the Birkhoff normal 
forms). Hence, the Birkhoff invariants should growth much more rapidly 
then the coefficient of the initial Hamiltonian and the set HB- m n A4 m is 
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sent to a much bigger domain by the mapping BNT m . Consequently, we 
cannot compute the measure /U< m of an image set in BMT m {HB- m n Mm) 
by the change of variables formula. 

In order to overcome this problem, we rescale the normalized variables in 
C 2n on the ball of radius s m around the origin where the Birkhoff transfor- 
mation $ m at order m is defined. Denoting the initial variables by z 6 C 2n , 
we consider z' € C 2n defined by z = & m (s m z') which correspond to a con- 
formal symplectic transformation giving the Hamiltonian : 

Km(z') = -^-H o $ m {s m z') or K m (x', y') = \-H o <S> m (s m x', s m y'). (5) 
s m s m 

In these rescaled variables, the integrable part k m (z') = -^-h m (s m I(x', y')) 

belongs to HB- m . 

Our initial issue is settled but the price to pay is a huge deformation 
of the measure fi< m since the Jacobian of the transformation BMT m which 

maps to the rescaled Birkhoff invariants (denoted ( Qi, • • • , Q m ) in the se- 



quel) is now s m where D(m) growth like dimension of the space of poly- 
nomial M m [X] of degree minn variables hence D(m) is of the order of m n . 

Let Xe™^ be a measurable set given by Yomdin quantitative Morse Sard 
theory of polynomial integrable Hamiltonians in HB- m which include the 
set of critical Hamiltonians at order m where Nekhoroshev theory cannot 
be applied in the rescaled variables to prove superexponential estimates of 
stability. 

We must bound the measure of the set of Hamiltonian in HB- m n M. m 
which are sent in by BNF m , hence : 

[l~(m)(BN'T m (P 1 , . . ■,P m ))dn< m (P) = ^ri) P-x(m)(Qi, ■ ■ ■ ,Qm)dfj,< m (Q) 

A s m J 

Mm BAfTm(Mm) 

and this is bounded by — — . 

With the growth of D(m), the divisor s m D ' m ' decreases in a superexpo- 
nential way with respect to the order of normalization m. At this step, the 
set introduced by Bounemoura is not big enough to ensure that its preimage 
in the construction of the Birkhoff normal forms is a generic set in % s . More 
specifically, we cannot afford a geometric decrease of the volume the critical 
parameters as in the convex case or in Bounemoura construction with the 
setxi m) . 

Here, we need a superexponential decrease of the volume of xi m \ 
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In order to settle this issue, we will not use Nekhoroshev theory on the 
whole domain where the Birkhoff normalization at order m is defined (hence 
on the ball of radius s m ) but on a much smaller domain which is the ball 
B rm of radius r m around the origin in the rescaled variables (hence r m s m in 
the normalized variables). 

After a normalization at order m £ N* and our rescaling, the size of the 
perturbation is of the order of (m!) 1+T (s m r m ) m ~ exp(— (1 + r)m)r™ (the 
prefactor comes from the growth of the coefficients). Hence the size of the 
perturbation is dominated by r™ since the expansion of the perturbation 
starts at order m (this is another important point), with the proposition 3.1 

and the previous reasonings we have ji< m of the order of (r m ) ma . 

Now we come back in our infinite-dimensional setting and, for all m £ N*, 
we consider the set X^^C M[[X]] such that G X^™ 1 ^ if and only if 

Proj Rm[x] (M is in (BNF m y\x { r ] ) or BNT m (vio] Wm[x] (h^eX^ . 
For all m <G N*, the product |~| is a probability measure and we 

k>m+l 

have n = n< m x fi k , hence : \i (xi m) ) = //< m \{BMT m ^ (xi m) ) 

k>m+l ^ 

With the previous reasonings for all m £ N*, we have the upper bounds : 



(x(«»)< 



D(m) D(m) 



Finally, the set of critical Hamiltonian C where Nekhoroshev theory 
cannot be applied at some order m 6 N* to prove superexponential results 
of stability is included in X^™ 1 - 1 which is of small measure if the sum 

raGN* 

(„ \ma 
D(m) 

The second radius r m is a /ree parameter which has to be chosen small 
enough to ensure the convergence of — Dtm) ' 

mtzfq* S m 

This imposes to choose a radius r m which is exponentially small with 
respect to the order of normalization m, here m is of the order of log(r m ) a 
where a > depends only of n. We cannot recover an exponentially small 
remainder with m of the order of l/r m as in the convex case but only a 
remainder which is smaller than any polynomially small remainder with re- 
spect to the radius. This is an important point which comes exclusively 
from the expected divergence of the Birkhoff normal form (actually, as pre- 
viously said, we know that this normal form is always obtained as a Gevrey 
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divergent expansion in the Diophantine case but we don't know if it diverges 
or converges generically) . Another misleading fact is that on the ball of ra- 
dius s m r m , we could obtain directly a Birkhoff normal form which integrate 
the system up to a super exponentially small remainder with respect to m. 
But in this case, the Birkhoff normal form would give only a remainder 
which is exponentially small with respect to the radius while here we start 
with a perturbation which is smaller than any polynomially small remainder 
with respect to the radius before going to the exponential with Nekhoroshev 
estimates. 

Summarizing we have proved that there exists positive constants a, Ci, 
C2, C3 and a sequence of annuli of decreasing radii (r m ) mS N* such that, for 
an arbitrary function H € Ji s and for almost all h £ HB^ with respect 
to the probability p, every solution of the Hamiltonian system governed by 
H+h with an initial condition r m+ \ < ||X(0)|| < r m in the rescaled variables 
satisfies : 

-C 3 \ln(r m )\ a 



Tit) - 1(0) <Cir m for\t\< exp C 2 r, 



rn 



In order to obtain uniform estimates with respect to the initial radius 
p, as in our theorem 1.4, it only remains to slightly lower the coefficients 
Ci, C2, C3 and the exponent a. 



4 Comments and prospects : 

A natural extension of this work would be for the case of a Lagrangian 
invariant torus for an Hamiltonian system. In this setting, we consider a 
system governed by the Hamiltonian u.I + /(I, <p)I 2 where / is analytic in 
action-angle variables (I, ip) 6 M. n x T n . The only missing ingredient to make 
the previous reasonings is a quantitative construction of the Birkhoff normal 
form to a given order in action. Up to the author knowledge, the available 
constructions ( |Fas90j , |P W94] . }Pos93j ) in the literature yield an integrable 
Birkhoff normal form with a perturbation which is small but with terms of all 
order in action including low order terms. Here, we need a Birkhoff normal 
form of the kind h m (I) + f m (I, *p)I m for m G N* to follow the scheme of our 
proof. This is a purely technical issue and there is certainly no obstruction 
to overcome this problem. 

The second extension would be for the case of an Hamiltonian with a 
lower regularity, that is for Gevrey or differentiable Hamiltonians. With the 
classical papers on Birkhoff normal forms in the differentiable case or the 
work of Popov-Mitev [MP10J in the Gevrey case together with the article 
of Bounemoura [Boullbj on generic Nekhoroshev estimates for Gevrey or 
differentiable Hamiltonian we have all the required technical tools. 
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A final interesting extension would be for the case of invariant tori pro- 
vided by KAM theory. A model statement would be : for an analytic inte- 
grable hamiltonian which satisfy Kolmogorov (or Russman) non- degeneracy 
condition ([AKN97J), KAM invariant tori are super exponentially stable for a 
small enough generic perturbation. The papers of Popov (|PopOO|, |Pop04|) 
on the existence of a global Birkhoff normal form for all the tori given by 
KAM theory should be useful in this setting. 
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